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I. INTRODUCTION 



After more than three decades of study, the number of proposed particle physics solutions 
to the solar neutrino problem is still increasing with time. The currently viable solutions 
to the available set of experimental data include two, three, and four neutrino oscillation 
scenarios (with vacuum and MSW oscillations among active as well as sterile neutrinos), 
neutrino decay, violation of Lorentz invariance, violation of the weak equivalence principle, 
and magnetic-moment transitions. Even for the simplest case of two neutrino oscillations, 
there are several isolated regions in neutrino parameter space that are consistent with all of 
the published data by the chlorine, Kamiokande, SuperKamiokande, GALLEX, and SAGE 
experiments. 

The existing solar neutrino data provide at most (2-3) a indications favoring specific 
solutions. Moreover, the predicted sizes of those neutrino conversion effects that do not 
depend upon the standard solar model and that can be measured well in the Sudbury 
Neutrino Observatory (SNO) [|I|], are typically small: from a few per cent to about ten per 
cent [H. Exceptions include the day-night asymmetry (for limited values of the oscillation 
parameters) and the double ratio of the neutral- to charged- current event rates. We will have 
to be lucky for the oscillation effects to be realized near their maximal possible values. In the 
largest part of the Am^ — sin^ 29 parameter space, currently acceptable neutrino oscillation 
scenarios predict that most of the new physics effects for SNO will typically be less than 
3 or 4 0" different from the no-oscillation predictions. And, as previous experience teaches 
us. Nature seems to prefer toying with us by providing ambiguous hints. The existence of 
systematic effects at the several percent level further increases the difficulty of identifying a 
unique solution. 

In this paper, we show how the predictions for some solar neutrino observables are 
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correlated, or why they are uncorrelated, in the context of specific solutions of the solar 
neutrino problems. We demonstrate that the signatures of a given neutrino solution include 
not only the values of specific observables, but also the correlations among the observables. A 
study of the correlations (and, where relevant, the lack of correlations) between the different 
predicted values of neutrino observables can be used to increase our understanding of the 
physical processes that are occurring.0 

In addition, there are excluded regions that we call "zones of avoidance." None of the 
currently favored oscillation solutions predict that the measured neutrino observables will 
lie within these regions of avoided parameter space. We stress the diagnostic value of testing 
predictions that new measurable quantities lie outside these current zones of avoidance. 

The main point of this paper is that studying the predicted correlations and zones of 
avoidance among solar neutrino observables can add discriminatory power to solar neutrino 
experiments. Although we illustrate the methodology using SNO variables and a partic- 
ular set of allowed neutrino oscillation parameters, the diagnostic value of the predicted 
correlations and zones of avoidance are more general than the particulars of our illustrative 
calculations. 

A. Previous discussions of correlations 

Correlations of solar neutrino observables have been discussed in several previous inves- 
tigations. Perhaps the first such discussion pointed out (Ref. 0) the lack of consistency, if 
no new physics were involved, between the total rates observed in the chlorine and in the 
Kamiokande experiments. This inconsistency was used as an argument to exclude astro- 
physical solutions of the solar neutrino problem. Kwong and Rosen (Ref. 0) analyzed the 
relations between event rates measured with SuperKamiokande and with the Sudbury Neu- 



*Also, numerical codes for calculating neutrino oscillation processes can be tested by requiring 
that they yield correlations predicted by analytic arguments given in this paper. 
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trino Observatory (SNO) for various MSW solutions. Even more closely related to what we 
discuss in the present paper, Folgi, Lisi, and Montanino ||^ mapped the Large Mixing Angle 
(LMA) MSW and the Small Mixing Angle (SMA) MSW solution regions in the Am^-sin^ 29 
plane onto the plane of flux independent observables: the day-night asymmetry, and the shift 
of the first moment of the electron spectrum measured with Super Kamiokande and SNO. 
The goal of the Bari group was to show how correlations could be used to help distinguish 
between the SMA and the LMA solutions. The correlations of a spectrum distortion and 
the day-night asymmetry for SMA and MSW Sterile solutions were discussed in Ref. 

As discussed in Ref. 0, strong correlation exists between the day-night asymmetry and 
the seasonal variations of signals for the MSW solutions: both effects originate from the 
same earth regeneration effect. For the SMA solution, a strong correlation exists between 
the total day-night asymmetry and the event rate in the "core bin" (the night bin in which 
neutrinos cross the core of the earth) [|10|. For vacuum oscillation solutions (VAC) it has 



been pointed out in Ref. [|ri|] that there is a strong correlation of spectrum distortion and 
seasonal variations. 

In Ref. [0 , we described the goal of analyzing simultaneously all of the SNO observables, 
measured and upper limits, in order to best constrain the allowed neutrino solutions. As a 
first step in that direction, we considered some pairs of measurable quantities but did not 
calculate the correlations between the predictions in the planes formed by the observables. 

In this paper, we illustrate the power of studying the correlations between different solar 
neutrino observables by evaluating the correlations between measurable quantities in the 
SNO experiment. We elucidate the physical basis for the correlations with the aid of simple 
analytic approximations. 

B. Correlated SNO Observables 

We consider here the correlations between the following quantities that are measurable 
in the SNO experiment. 
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The total reduced rate of the charged- current events above a specified threshold: 

[CC] - ^ . (1) 

CC 

Here A^^cc is the observed number of events from the CC (neutrino capture reaction 
by deuterium) and A^cc^ number expected on the basis of the BP98 standard 

solar model and no new particle physics beyond what is predicted by the standard 
electroweak model. The predictions for [CC] implied by the six currently allowed two- 
neutrino oscillation solutions have been calculated in and are also discussed in 



In what follows, we consider first the correlations of different experimental quantities 
with [CC] since the CC rate is the easiest quantity for the SNO collaboration to 



measure. 



The day-night asymmetry of the charged- current events 14 : 



CC _ N-D 

^N-D = . (2) 

Here N and D are the rates of the events observed during the night-time (N) and 
the day-time (D) averaged over the year (and corrected for the seasonally changing 
distance between the sun and the earth). The contours of constant v4^^j-, have been 
calculated in Ref. [0. In what follows, we will use the notation v4n_e) to denote the 
charged- current day-night effect and will use the more cumbersome notation, y4^^j-), 
only when there is a chance of confusion with the day-night effect measured from 
neutrino-electron scattering in SuperKamiokande, ^n-d- 

The relative shift of the first moment of the electron energy spectrum from its non- 
oscillation value : 

6T^^. (3) 
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Here T and Tq are the first moments of tfie recoil electron energy distribution calculated 



with and without neutrino oscillations. The shift has been defined in [0; we calculated 
6T in Ref. for the currently allowed set of neutrino parameters. The distortion is 
expected to be smooth for all solutions except for vacuum oscillations with large Am^, 
so that the first moment characterizes the distortion of the recoil energy spectrum 
rather well. 

The double ratio of the reduced neutral-current rate (neutrino disintegration of deu- 
terium), A^NC; to the reduced charged-current event rate: 



[NC] _ iVNc/iVcc 



[CC] {N^c/Ncc] 



SSM 



(4) 



We will also discuss the ratio of the reduced rates of neutrino-electron scattering [ES] 
and charged current events [CC]: [ES]/[CC], 



[ES] _ Nes/Ncc 



[CC] - {N^s/Nccr^^ 



(5) 



Here [ES] = A^es/^Is^' where A^es is the number of observed u — e scattering events 
and A^lg^ is the number of predicted events according to the SSM. 

Additional SNO observables are discussed in [E|. In particular, the seasonal variations 



may be significant for both vacuum |17| and MSW |0,^8[ solutions 



C. Outline of this paper 

In Sec. 0, we describe our method. In the next three sections, we study correlations 
related to the charged-current (CC) events observable with SNO: An-d— [CC] in Sec. fT|, 
[CC] - 6T in Sec. fV|, and An-d — ST in Sec. |V|. We discuss in Sec. |V3 correlations in 
the plane of [NC]/[CC] and An-d and in Sec. [VT] the correlations of [NC]/[CC] and 5T. 
In Sec. |V111| , we summarize our results. In the Appendix, we describe the dependence on 
oscillation parameters of each of the observables discussed in the main text. Simple analytic 
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expressions for these dependences are presented in the Appendix; these analytic expressions 
are useful for interpreting the results of detailed numerical calculations. 



D. How should this paper be read? 

The most efficient way to read this paper is to first obtain an overview of what is ac- 
complished and then to descend into the details. We recommend that the reader begin by 
looking at Fig. ^ to Fig. |ll], which show the correlations and the zones of avoidance in planes 



constructed from quantities that are measurable by SNO. Then we suggest that the reader 
turn immediately to Sec. |VIII| , where we provide a succinct summary of our principal results 



and conclusions. Only after having acquired this overview, do we recommend returning to 
Sec. H in order to read about the details. 

II. METHOD 

For specificity, we consider correlations of observables in the SNO experiment for the 
two-fiavor neutrino solutions of the solar neutrino problem. Each solution is characterized 
by the two oscillation parameters, /S.rn? and sin^ 26. We use the techniques described in 



Ref. to determine the allowed regions for the oscillation parameters. The input data 
used here include the total rates in the Homestake, Sage, Gallex, and SuperKamiokande ex- 
periments, as well as the electron recoil energy spectrum and the Day-Night effect measured 
by SuperKamiokande in 825 days of data taking. 

Fig. HI shows the acceptable regions of the solutions in the plane of /S.rn? — sin^ 26 as 



determined in Ref. |jT3[ . For our study of correlations, as exhibited in Fig. 0-Fig. |TT 



we use 



the 99% CL. solutions shown in Fig. |I[ 

We stress that the particular topography of the predicted correlations and the zones of 
avoidance depend upon the set of experimental data that are used in finding the allowed 
regions and the Confidence Limit (CL) that is adopted. We have adopted the input data 
and the CL specified above. Of course, the correlations and the zones of avoidance will 
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FIG. 1. Global oscillation solutions. The input data include the total rates in the Homestake, 
Sage, Gallex, and SuperKamiokande experiments, as well as the electron recoil energy spectrum 
and the Day-Night effect measured by SuperKamiokande in 825 days of data taking. Figure |l|a 
shows the global solutions for the allowed MSW oscillation regions, known, respectively, as the 
SMA, LMA, and LOW solutions [Q. Figure |l|b shows the global solution for the allowed vacuum 
oscillation regions. The CL contours correspond, for both panels, to = Xmin + 4.61(9.21), 
representing 90% ( 99% CL) relative to each of the best-fit solutions which are marked by dark 
circles. 

evolve as more experimental data become available. Hopefully, the predicted correlations 
will become stronger and the zones of avoidance much larger. 
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We perform both numerical and semi-analytical studies of correlations. The analytic 
work provides, in many cases, a simple interpretation of the numerical results. 

For SNO, two solar neutrino fluxes are relevant: the flux and the hep flux. We 



characterize these two neutrino fluxes by the dimensionless parameters /s and fhep pO| , |2T 
which are the fluxes in units of the BP98 Standard Solar Model fluxes 



A. Mapping from neutrino space to observable space 

Predicted values of SNO observables, X (e.g., [CC], An-d, [NC]/[CC]), are functions 
of two oscillation parameters and two flux parameters: 

X = X{^m\^iTi^2ejBjhep). (6) 



Following the same procedure as in Ref. [0, we determine fs and fhep for each pair of 
values of the oscillation parameters, Am^ and sin^ 26, by fltting to the total rate and the 
recoil electron energy spectrum of SuperKamiokande ||22||: 



fB = fB{^m\ sin^ 29), hep = hep{^m\ sin^ 29). (7) 

After this determination, the SNO observables are functions of two oscillation parameters 
only: 

X = X(Am^ ?,in^29). (8) 

The correlations depend on the form of the functions, Eq. (H); the functions are different for 
each solution of the solar neutrino problem. We give in the Appendix simple parametriza- 
tions of the dependences for various oscillation solutions. 

In what follows, we flnd regions in planes of (X, Y) observables allowed by the data from 
all existing solar neutrino experiments. Formally, this is equivalent to mapping the solution 
regions in the /S.rn? — sin^ 29 plane onto the plane of observables oi X - Y . For each point 
Am^ — sin^ 29 of the solution regions, we calculate values of X and Y . The mapping is given 
by Eq. (I). 
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If the region of an oscillation solution in the Am^ — sin^ 29 plane is projected onto a hne 
segment or onto a narrow strip in the X and Y space, then we say that there is a strong 
correlation of the observables X and Y for the specified solution. In some cases, there is a 
strong correlation only in part of a given solution region. 

There are various ways one might quantify the degree of correlation. The most appropri- 
ate way for our purpose (enhancing the identification power of the analysis) is the following. 
Let us denote by Sxy the area of the region in the X — Y plane to which a given solution 
region is projected. Let /S.X and Ay be the intervals of the observables X and Y in which 
these observables can vary within a given solution region if we consider the variables as 
independent. The product AX x Ay is the area of the mapped region if X and Y are 
uncorrelated. The degree of the correlation of X and Y can be characterized by the ratio: 

'^^^ = AXxAy • 

If the correlation parameter kxy << 1, we will say that a strong correlation of the X and Y 
observables exists. In this case, the allowed X and Y parameter space is small and "zones of 
avoidance" dominate. For strong correlations, a combined study of the observables X and 
y will enhance the identification power of the analysis. If kxy ~ 1, there is no correlation 
and no advantage to a combined study of X and Y . 

B. Analytic approximations 

The accurate prediction of solar neutrino observables requires multiple integrations over 
energy-dependent survival probabilities and neutrino interaction cross sections, and also 
over the energy resolution and the efficiency of detection. In spite of the complicated nature 
of the accurate calculations, simple and useful analytic results can often be found. The 
analytic expressions generally contain a small number of parameters that can be determined 
using the detailed numerical results. In developing analytic approximations, we proceed as 
described below. 
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First, we determine the functional dependence of the observables on the neutrino oscil- 
lation parameters, primarily through the dependence of the survival probability, P, on the 
oscillation parameters. Thus 

X{Am\sm^2e) - X{P{Am\sm^2e)) , (10) 

where the parameter P represents the survival probability after a suitable average over the 
energy. Therefore the first step is to find expressions for observables in terms of P. 

Second, the expressions for the survival probability P can usually be simplified in the 
restricted regions of oscillation parameters that apply to specific allowed solutions. Also, 
averaging over relatively small intervals of energies (smoothing the dependences) often leads 
to further simplification. 

Third, in the analytic expressions for observables, the average neutrino energy should be 
taken as a fitting parameter, which is determined by comparison of the analytic expression 
with the detailed results of numerical calculations. The energy parameter should be fitted 
separately for different solutions and for different observables. Moreover, if a given observ- 
able is described by several terms with different dependences on energy, the characteristic 
energy in each term should be considered as an independent parameter. This procedure 
will usually give the correct parametrization provided that there are no particularly strong 
energy dependences. The approximation generally works well if the fractional change of 
the survival probability over the effective range of integration is reasonably small. This 
condition is usually satisfied for most SNO observables. 

In some cases, e.g., when observables depend on the same combination of oscillation pa- 
rameters, the exact results for correlations can be obtained without performing complicated 
integrations over energies and over instrumental characteristics. 

In summary, we find dependences of observables on oscillation parameters in terms of 
simple functions with a few fitting parameters that are determined by exact numerical 
calculations. The fitting parameters represent the complicated results of integrations over 
energies and over instrumental characteristics. 
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C. Survival probabilities, observables, and correlations 



We find in this subsection the dependence of different neutrino observables on the average 
survival probability. This is the first step in the process of deriving analytic expressions 
for correlations, which was outlined in the previous subsection. We shall see that some 
correlations appear clearly even when only survival probabilities are considered. In the 
Appendix, we derive expressions for the survival probabilities and show how these expressions 
can be used to predict correlations among neutrino observables. 



1. Charged current in SNO and neutrino- electron scattering in Super Kamiokande 
The reduced CC-event rate in the SNO detector can be written as 

[CC] = PsNO ■ fB , (11) 

where Psno = -Psno(^^^5 sin^ 2^, is the effective survival probability for CC events in 
SNO experiment. The hep neutrinos do not contribute significantly to the total rate for an 
energy threshold in the likely range of 5 to 8 MeV and therefore hep neutrinos can usually 
be neglected. We determine the flux parameter fs from the reduced total rate of events in 
the SuperKamiokande neutrino-electron scattering experiment: 

[ES]sK = ATsK/iVll^, (12) 



where A^'sk and A*"!^^ are the observed and the SSM (BP98, see Ref. [0) predicted event 
rates, respectively. In the case of oscillations into active neutrinos we find 

f ^ [ES]sK 

-fSK + U — -fSKj'^ 



where r ^ 0.16 |2^ is the ratio of the z/^ — e to the z/e — e cross-sections, and Psk 
Psk(^^^; sin^ 26*, E^^) is the effective survival probability for the SK experiment. 
Thus, we find from Eq. (0) and from Eq. (|13D 



^ J ^ (l-r)(PsK/PsNo)+r/PsNO ' ^ ^ 
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where the ratio Psk/-Psno depends in general on the oscillation parameters but is often of 
the order of unity. Equation (|T^) simplifies considerably if we make the approximation (valid 
for example if the SK and SNO energy thresholds are chosen near the plausible values of 
~ 6.5 MeV and E^^^ ~ 5 MeV, see Ref. 0) that PgK ~ i^sNO = P- In this special 



circumstance, 

[CC] ^ ^^^^^^ . (15) 
^ ^ 1-r + r/P ^ ^ 

Equation (|I5|) is generally valid for the LMA and LOW solutions 0, for which the survival 
probability depends only weakly on the energy in the energy range of interest. 
In the case of conversion to sterile neutrinos (r ^ 0), we have 



[CC]sterile = ^.^^ ^ m,^- (16) 
-rsK 



For PsK ~ PsNO, the rate [CC] is approximately equal to [ES]gj^. 

Deviations from the equality Psk = -Psno can be caused by a strong energy dependence 
of the survival probability, by differences in the energy dependences of the neutrino cross- 
sections, by difference of energy thresholds, and by differences in instrumental responses. 

2. Shift of the first moment of the CC spectrum 

The fractional shift of the first moment of the recoil electron energy spectrum [see Eq. (^] 
is easily shown (for a negligible hep flux) to be proportional to the derivative of the survival 
probability: 

""PIE^ ^'^^ 
where P and dP/dE are suitable spectrum averages and ii^ is a characteristic energy. 

3. Day-Night asymmetry 
The day-night asymmetry can be estimated from the expression 
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where Pn and Pd are the day and the night survival probabihties averaged over the year 
after removal of the geometrical factor R~^. We showed previously in Ref. that the day- 
night asymmetry in the SNO CC-events, ^n-d? the asymmetry in the neutrino-electron 
scattering observed by SuperKamiokande and SNO, ^n-d; related by the approximate 
equation 



aCC _ 4ES 



Equation ( [T9| ) shows that the CC day-night asymmetry that will be measured in the SNO 
experiment is predicted to be larger than the neutrino-electron scattering asymmetry mea- 
sured by SuperKamiokande. For typical values of the average survival probability in the 
range P = 0.3 — 0.5, the enhancement factor in brackets in Eq. (|19D is between 1.4 and 1.6. 
The prediction given in Eq. (|^) can be tested also by using SNO data alone since both 
A^'^Y) ^iid A^^Y) measurable by SNO. 

Combining Eqs. (|l^) and ([T9|) , we obtain a relation between the day-night asymmetry 
^N-D the reduced CC rate [CC] : 

[ES] 

SK /o^^ 

^N-D - [CC](l-r) ' ^^^^ 

where ^n-d depends mainly on the Am^ for the LMA and LOW solution regions (see 
Ref. PI): ^N-D = ^n-d(^"^^)- This relation holds pointwise, i.e., for a particular choice 
of Am^ and sin^ 29. Most of the range in the ^n-d [CC] plane is due to the allowed 
range in Am^ and sin^ 29, which washes out the pointwise dependence of Eq. (|20| ) because 
in the LMA and LOW solution regions ^n-d depends primarily on Am^ and [CC] primarily 
depends upon sin^26' (cf. Fig. |^ and Fig. ^ ). 

4. The double ratio [NC]/[CC] 



The double ratio [NC]/[CC] is equal to the inverse of the appropriately-averaged survival 
probability for the active neutrino case: 
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[NC]/[CC] = i . (21) 



Both [NC]/[CC] and [CC] are determined by P [see Eq. (|T|)]. Inserting Eq. (|21|) into 
Eq. (P^D, we obtain 

f^^] l-r|^r[NC]/[CC] ' ^^^^ 
which imphes that [NC]/[CC] and [CC] are strongly correlated in our approach (in which 
fs is fixed by the measured SuperKamiokande rate). As a consequence of Eq. (|22D , the 
correlation plots are similar for [NC]/[CC] and [CC] when combined with other observables. 
A strong correlation exists also between the double ratios [ES]/[CC] and [NC]/[CC], both 
of which will be measured by SNO: 

M-i_, + ,M (23) 

[CC] " ^ ' + ' [CC] • ^^^> 

For sterile neutrinos, 

[NC]/[CC]3,,,,, = ? ^ ^ ' (24) 

where P' is the average survival probability for the NC event sample. Since the thresholds 
for NC events (2.2 MeV) and for CC events (expected to be greater than 5 MeV) are 
different, the ratio of the average survival probabilities, P' / P, is in general different from one. 
However, for both NC and CC events the cross section increases with neutrino energy and 
most of the events that are observed correspond to neutrinos with relatively high energies. 
For these higher energy neutrinos, the survival probability depends rather weakly on energy 
(see Fig. 1 of Ref. [0]). So, while P'/P is not identical to one it is in general quite close to 
one for practical cases. 

Combining Eq. (|T^ and Eq. (^), we find 

'1 — r' 



[NC]/[CC] 



4CC _ 4ES 

^N-D ^N~D 
^N-D 



(25) 



Equation ( pSj ) is an example of a correlation between three observables. The equality given in 
Eq. (pSD does not depend explicitly on the oscillation parameters and holds approximately 
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for all three MSW active neutrino solutions. The principal inaccuracy introduced in the 
derivation of Eq. p5| ) is caused by the fact that the average survival probability, P, that 
appears in Eq. (|I9|) is not exactly equal to the average survival probability that appears in 
Eq. (^. 

If we neglect the dependence of the measured quantities upon energy threshold and 
upon the energy dependence of the rates, then the neutrino observables depend only on two 
parameters, /S.rn? and sin^ 20 . Therefore any three observables X, Y, Z must be correlated, 
except for special cases in which one or more of the observables do not depend on the 
oscillation parameters. Indeed, expressing the oscillation parameters in terms of the two 
observables, say X and Y, we can get the relation Z = Z(Am^, sin^ 20) = Z{X,Y). The 
experimental study of the validity of Eq. (]25|), and other similar "triple" relations, will 
provide important tests of the consistency of the oscillation solutions and the experimental 
results. Deviations from the "triple" relations that could not be explained by expected 
energy dependences of the experimental quantities, or by differences in the average values 
of P for the various measurables, would indicate either the participation of more than two 
neutrinos in solar neutrino oscillations or a lack of consistency of the experimental results. 

5. What's next? 

The functional dependence of the survival probability on the oscillation parameters de- 
pends on which particular solution of the solar neutrino problems is chosen. In the Appendix, 
we give the function dependences for different currently-favored oscillation scenarios. Using 
the expressions for P given in the Appendix and the relations presented in Eq. ([151) , Eq. (p!7D, 
Eq. (|18D, and Eq. (pID, we derive the dependences of the SNO observables on the neutrino 
oscillation parameters. 

In the next three sections, we present maps of neutrino oscillation solution regions onto 
planes constructed from different pairs of SNO observables. We discuss results for the 
following currently-favored two-neutrino solutions which explain all of the available solar 
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neutrino data: Large Mixing Angle (LMA) MSW solution, Small Mixing Angle (SMA) MSW 
solution, Low Am^ (LOW) solution, and MSW Sterile solution based on small mixing angle 
MSW conversion to sterile neutrinos. There are several disconnected regions ("islands") of 
the Vacuum Oscillation solutions. We will divide them into two groups: vacuum oscillation 
solutions with small Am? {Am? < 10^^'^ eV^), VACs, and several islands with large Am^ 
(Am^ > 10"^° eV^) solutions, VACl. For VACs, there are currently four allowed islands in 
the Am^ and sin^ 26 plane. The VACs solutions correspond to four almost fixed values of 
Am^ and varying sin^ 26. 



III. CC-RATE VERSUS DAY-NIGHT ASYMMETRY 

Figures and ^ show maps of the currently- allowed regions in the Am^— sin^ 26 plane onto 
the [CC] — An-d plane for the electron energy thresholds of 5 MeV and 8 MeV, respectively. 
In Fig. 0, we show a simulated data point near the best-fit value for the LMA solution. 
The estimated la error bar for the CC measurement is taken from Table II of Ref . [0] . For 
An-d, "we assume for purposes of illustration a ±0.03 uncertainty in the absolute value as 
a la error, which is comparable with the accuracy that has been achieved after three years 



with the SuperKamiokande detector 22]. 



A. Discriminating among solutions 

There are four regions in Fig. ^ and Fig. ^ for which (after taking account of the likely 
measurement uncertainties and the overlap of the predicted values of the observables) im- 
portant scientific inferences will be possible if the measured values of [CC] and An-d fall 
within the designated areas. 1) If measurements show that A^^d > 0.2, then that will be 
a strong indication in favor of the LMA solution. 2) If the measurements only show that 
^N-D > 0.1, that by itself will be sufficient to disfavor the vacuum and MSW Sterile solu- 
tions. 3) If instead Ad-n ~ and [CC] is consistent with 0.48, then that would strongly 
favor the MSW Sterile solution. 4) If the measured values lie within the 'zone of avoidance' 
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FIG. 2. The allowed regions for the day- night asymmetry in the charged current event rate, 
j4n-D5 versus the reduced charged current rate, [CC], for an electron energy threshold 5 MeV. The 
figure shows the currently-allowed regions predicted by two-neutrino solutions ||l^ that describe all 
the available solar neutrino data: LMA (encircled by a solid line), SMA (dashed line), LOW (dotted 
line), VACs (black points), and VACl (grey points). The best-fit points for each solution are 
indicated by a small black circles within the allowed region. The prediction for the no-oscillations 
case is indicated by a triangle. The cross near the best-fit point of the LMA solution is a simulated 
measurement with estimated la error bars. 

of An-d > 0.02 and [CC] > 0.4, then none of the currently acceptable oscillation solutions 
will be favored. 

If An-d > 0.1, then the inferences from cases 1) and 2) above can be tested by measuring 
6T. Figure § and Fig. H, as well as Fig. || and Fig. ^, show that all the currently allowed 
oscillation solutions predict 6T < 0.01 if An-d > 0.15. 

If the day-night asymmetry lies in the broad range < An-d < 0.15, it will be difficult 
to disentangle the LMA, LOW and SMA solutions. These three MSW solutions show a large 
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FIG. 3. The allowed regions for the day-night asymmetry, ^n~Di versus the reduced charged 
current rate, [CC], for an electron energy threshold 8 MeV. The meaning of the symbols is the 
same as for Fig. ||, except that the regions now refer to a recoil electron energy threshold of 8 MeV. 



overlap in their predictions for the [CC]-An-d plane (see Fig. |^ and Fig. |] and especially 
Table IX of Ref. 0). If the observed asymmetry is not too small, e.g., if A^-d > 0.1, then 
the SMA solution can be identified by the zenith angle dependence of the rate during the 
night. According to the SMA scenario, the rate should be strongly enhanced in the deepest 
night bin (for the core-crossing neutrino trajectories) In contrast, the LMA and LOW 
solutions predict rather flat zenith angle distributions. The LOW and the LMA solutions 
may be distinguishable through the observed dependence of the day-night asymmetry on 
the energy threshold. The asymmetry increases with threshold for the LMA solution and 
decreases with threshold for the LOW solution. For the LMA solution, the maximal possible 
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asymmetry becomes as large as 0.32 for = 8 MeV instead of 0.28 for E = 5 MeV (see 
Fig. and Fig. |^). For the LOW solution, the dependence upon the threshold energy is 
just the opposite; the predicted asymmetry decreases with increasing threshold energy. In 
particular, the LOW solution predicts that the maximal asymmetry decreases from 0.135 
to 0.11 as the threshold energy increases from 5 MeV to 8 MeV. The LOW solution may 
also be identified later by strong Day-Night variations of the beryllium line in BOREXINO 



experiment 



The charged-current event ratio is in some ways the simplest experimental quantity to 
measure with the Sudbury Neutrino Observatory. However, the most remarkable aspect of 
the above analysis of Fig. Q and Fig. |^ is that the potentially important inferences are almost 
entirely independent of the measured charged-current rate. This is because the estimated 
la uncertainty in the value of [CC] is about 6.7% (see Ref. 0) and is dominated by the 
theoretical uncertainty in the charged-current neutrino-absorption cross section. Unless a 
major improvement is made in the accuracy of the theoretical cross section calculation, the 
potential diagnostic value of the charged-current measurement will be severely compromised 
by the large uncertainty in the neutrino absorption cross section. 



B. Correlation phenomenology 

For the LMA solution, there is no significant correlation shown in Fig. ^ and Fig. ^; the 
correlation parameter, cf. Eq. @, fi;^_cc (LMA) ~ 1. For the largest area of the plane of 
oscillation parameters, the charged-current rate depends mainly on sin^ 26' [see Eq. ( |A8|) ], 
whereas An_£) depends strongly on Am^ [see Eq. ( [AlOD ]. There is a tendency for small 
values of [CC] to correspond to large values of An-d, since [according to Eq. (^)] for fixed 
Am^ the asymmetry is inversely proportional to [CC]. 

Also for the LOW solution, no significant correlation appears. The area occupied in the 
[CC]-y4N-D plane by the LOW solution is substantially smaller than for the LMA solution, 
which reflects the smaller allowed region of the LOW solution in the Am^ — sin^ 29 plane. 
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The SMA solution has the form of two beautiful, asymmetric petals connected at the 
point of zero asymmetry. This form can be understood from the expression for the asym- 
metry, Eq. ([A14|). The zero asymmetry contour is determined by the condition P = 1/2. 



Therefore, the contours of An-d = and of [CC] = 0.41, both of which correspond to 
P = 1/2, coincide. The contours are defined by the relation 

^ = Am^ ■ sin^ 29 = constant . (26) 

The correlation between An-d and [CC] appears in the region of small asymmetries and of 
large [CC]. The rate [CC] decreases with increase of An-d (Here /t^_C(7(SMA) ^ 1 ) [|. 

For vacuum solutions, there is a correlation between An_d and [CC] that is difficult to 
see on the scale of Fig. and Fig. ^ because the day-night asymmetry is small. The residual 
asymmetry, which is calculated after first removing the dependence of the total fiux, 
is not zero, but An-d < 2% is predicted for all the currently-favored vacuum oscillation 
solutions ||^. By inspecting the figures carefully, one can see that the asymmetry v4n_d 
increases with decreasing [CC]. 

The day-night effect for vacuum oscillations is determined by geometrical factors. In 
the northern hemisphere, the nights are longer in the winter when the earth is closer to the 
sun. The earth-sun distance affects the vacuum oscillation probability and therefore the CC 
event rate. Combining Eq. (15), Eq. (|A22|) , and Eq. ( A24 ), we find 



r 

where 



^N-D = -f^-llx/(Am^), (27) 



V [CC] 



/(Am^) ^ 2(AmVm^) cot(AmVm^) (28) 

is a function of Am^ only. Since within a given currently-allowed "island" in the 
Am^ — sin^ 26* plane , the variation of Am^ is small, we can consider /(Am^) ^ constant. 



''^A similar plot for correlation of the slope parameter and the asymmetry have been given in 
Ref. il. 
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Equation ( ]27|) explains the correlation between An~d and [CC] that exists in Fig. | and 
Fig. |. 



IV. CHARGED-CURRENT RATE VERSUS SHIFT OF FIRST MOMENT 

Figures ^ and |^ show, for electron energy thresholds of 5 MeV and 8 MeV, maps onto the 
[CC] — 6T plane of the currently-allowed regions in the Am^ — sin^ 29 plane. For illustrative 
purposes, Fig. | shows a simulated experimental point near the current best-fit predicted 
point for the LMA solution. The error bars are estimated la uncertainties from Table II of 
Ref. [3, with a 1.3% fractional uncertainty in the first moment and a 6.7% for the charged- 
current rate. 

A. Discriminating among solutions 

If the measured values of [CC] and 6T fall close to the current best-fit value of the 
LMA solution, then many of the currently-favored solutions will still be allowed if a 3a level 
of disagreement is permitted. The difficulty in uniquely identifying solutions is primarily 
caused by the estimated 3a uncertainties being comparable in many cases to the size of the 
predicted effects. 

There are some regions of the two-dimensional parameter space, [CC] — ST, that are 
relatively discriminatory. For example, in the region in which 6T > 0.04 and 0.3 < [CC] < 
0.4, only the VAC^ and the SMA solutions are represented. The most extreme values of the 
[CC] parameter, e.g., [CC] > 0.5 or [CC] < 0.3, would indicate, respectively, the MSW Sterile 
solution or the LMA solution. If either of these cases is suggested by the [CC] measurement, 
then a comparison of the predicted and measured 6T (Fig. ^) and An-d (Fig- 0) will be 
useful checks of the validity of the identification of the solution. 

Unique inferences will be possible (see Fig. ^ for extreme VAC5 solutions with a frac- 
tional shift 6T > 4.5% and for the MSW Sterile solution with [CC] greater than 0.48. The ex- 
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FIG. 4. The allowed regions for the shift of the first moment, 5T, versus the reduced charged 
current rate, [CC], for a recoil electron energy threshold of 5 MeV. The meaning of the symbols 
is the same as in Fig. ^, except that the regions now refer to ST and [CC]. The nearly horizontal 
segments that overlap with the LOW and the LMA solution regions correspond to VACs solutions 
with the largest values of Am^. 

treme VACs solution predicts a very small value for the day-night asymmetry, Ad-n < 0.01 
(see Fig. P). 

Two zones of avoidance appear in Fig. ^ for large [CC]. None of the currently favored 
oscillation solutions predict 6T < 0.01 or 6T > 0.04 for [CC] > 0.45. For an electron recoil 
energy threshold of 8 MeV, we find that 6T < 0.04 for all the currently favored oscillation 
solutions (see Fig. P). For smaller [CC], between 0.3 and 0.4, there is also a zone of avoidance 
for 6T > 0.01 and less than the values predicted by the SMA solution. 
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FIG. 5. The allowed regions for the shift of the first moment, 5T, versus the reduced charged 
current rate, [CC], for a recoil electron energy threshold of 8 MeV. The meaning of the symbols is 
the same as in Fig. ||, except that the allowed regions refer to 6T and [CC] and the threshold for 
the recoil electron energy is 8 MeV. 

B. Correlation phenomenology 



For the SMA solution, both [CC] and 6T are determined by a unique combination of 
neutrino variables, ^, defined by Eq. (p6|), so that (up to small earth matter effect corrections) 
the two measurables are strongly correlated. As 6T increases, the charged- current rate [CC] 
decreases (cf. Fig. ^). Using Eq. (|I4|) (with Psno ~ -Psk) and Eq. ( [A13| ) for the rate and 
Eq. ( |A16| ) for the shift of the first moment, we find for an electron threshold energy of 5 
MeV 



[CC] 



[ES] 



SK 



1 — r + r ■ e 



BST 



(29) 
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The numerical coefficient B in the exponent that occurs in Eq. (^) is determined by results 
of exact numerical calculations, B = 29.3. For an electron energy threshold of 8 MeV, one 
should use the general formula given in Eq. ( [T^ ) without making the approximation that 
-PsNO ~ -PsK- 

For the VACs solution, a strong correlation exists. Using Eq. (|T^) and Eqs. ( [A23D and 
(PI), we find 



= i (Hp - l) X f(Arn') . (30) 

where the function /(Am^), has been defined by Eq. (|28|) . The relation given in Eq. (|30D 
holds for each allowed island in neutrino parameter space, and, in the approximation of 
constant Am^, there is a strong correlation of the rate and the shift of the first moment 
for each of the three islands with low Am^. The first moment shift, 6T, increases as [CC] 
decreases. For the allowed island with the largest value of Am^ (which overlaps for a 5 MeV 
threshold with the LMA and LOW solutions), the shift of the first moment is close to zero 
and the correlations are weak. 

For the MSW Sterile neutrino solution, the rate [CC] is strongly restricted by the mea- 
sured value of [ES]gj^, whereas 6T varies over a significant range. 

The LMA solution does not predict a strong correlation, as can easily be seen from 
Eqs. ( |A8| ) and ( |A9D . The rate [CC] depends strongly on sin^ 29, whereas 6T depends strongly 
on Am^. The situation is similar for the LOW solution. 

At a higher threshold, 8 MeV (see Fig. |^), the shift of the first moment becomes smaller 
for all solutions. In particular, a significant part of the LMA region with negative 6T 
disappears. At the same time, for the VACl region the best-fit point shifts to larger 6T. In 
contrast, the spread of the [CC] rates increases especially for the SMA Sterile solution. For 
E^'^ = 5 MeV, PsNO ~ Psk and, as a consequence, [CC] is uniquely fixed by Rsk [Eq- (0)]- 
For i?*^ = 8 MeV, Psno/Psk differs from one and depends on oscillation parameters, which 
leads to the larger spread in [CC]. Also for E^^ = 8 MeV, the VACs region with the largest 
Am^ no longer overlaps with the LMA and the LOW solution regions. 
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FIG. 6. The allowed regions for the shift of the first moment, 5T, versus the day-night asymme- 
try, j4n-D) for an electron recoil electron energy threshold of 5 MeV. The meaning of the symbols 
is the same as in Fig. except that the allowed regions now refer to 6T and j4n-d- 

V. SHIFT OF FIRST MOMENT VERSUS DAY-NIGHT ASYMMETRY 

Figures ^ and ^ show maps of the currently-allowed regions in the Am^— sin^ 29 plane onto 
the 6T — A^_D plane for the electron energy thresholds of 5 MeV and 8 MeV, respectively. 
In Fig. 1^, the estimated 1.3% (la) error bar for the 6T measurement is taken from Table II 
of Ref. 0. For An-d, we assume for purposes of illustration a ±0.03 uncertainty in the 
absolute value as a la error, which is comparable with the accuracy that has been achieved 
after three years with the SuperKamiokande detector P^ . 
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FIG. 7. The allowed regions for the shift of the first moment, 5T, versus the day-night asymme- 
try, j4n-D) for an electron recoil electron energy threshold of 8 MeV. The meaning of the symbols is 
the same as in Fig. |2|, except that the allowed regions now refer to 6T and ^n~d and the threshold 
for the recoil electron energy is 8 MeV. 

A. Discriminating among solutions 



The only truly unique regions in the An-d — ST plane are the very large values of 
^N-D > 0.2, which would favor LMA, and the very large values of 6T > 0.06 (5 MeV 
threshold), which would favor VAC5. In both cases, the oscillation solution implies that the 
other measured parameter should be small, i.e., 6T should be small (according to LMA) if 
An-d is near its maximal value and An-d should be small if ST is near its maximal value. 

Both the LMA and the LOW solutions predict that the shift of the first moment is small 
for all allowed values of the day-night asymmetry, i.e., —0.02 < 6T < 0.015. And, of course, 
the day-night asymmetry is predicted to be small for all the vacuum solutions and the MSW 
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Sterile solution, i.e., I^n-dI < 0.02 for all allowed values of 5T. The imposition of these 
cross checks can be used to test the validity of currently-allowed oscillation solutions. 

Taking into account the estimated uncertainties in the measurements, the most populated 
region in the An-d — plane contains multiple currently-allowed solutions. 

It is easy to find zones of avoidance in Fig. ^ and Fig. |^. For a 5 MeV electron recoil 
energy threshold, there are no predicted solutions with An-d > 0.02 and 5T > 0.045 nor are 
there any predicted solutions with An-d > 0.1 and 6T < 0. For an 8 MeV energy threshold, 
there are no predicted solutions with An-d > 0.02 and 6T > 0.01. 



B. Correlation phenomenology 

In the small Am^ limit (large day-night asymmetry), both the LMA and the LOW 
solutions predict an approximately linear relation between the day-night asymmetry and 
the fractional shift in the first moment: 

6T = fc^N-D • (31) 

For the LMA solution, /clma = 0.014 and /^low = —0.03 for a 5 MeV threshold. These results 
can be obtained from Eq. ( [A9| ) and Eq. ( |A10| ) for the LMA solution and from Eq. (|A20|) 



and Eq. ( [A2 1|) for the LOW solution. This weak correlation exists because in the region of 



small Am^ of the LMA solution and large Am^ of the LOW solution both the day-night 
asymmetry and the shift of the first moment are induced by the earth matter effect. 

The allowed regions for the SMA and the MSW Sterile solutions both have the form of 
two petals connected at the point An-d = 0. 

For vacuum solutions, with the accuracy that is apparent in Fig. ^ and Fig. |^, the shift 
of the first moment does not seem to depend significantly upon the day-night asymmetry. 
However, there is a linear correlation, 

Ad-n oc ST, (32) 
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FIG. 8. The allowed regions for the [NC]/[CC] double ratio versus the day-night asymmetry, 
An-D) for a recoil electron energy threshold of 5 MeV. The meaning of the symbols is the same as 
in Fig. ^, except that the allowed regions now refer to [NC]/[CC] and ^n-d- 

where the coefficient of proportionahty is sufficiently small that the variation of about 
zero is not easily visible on the scale shown in the figures. The correlation arises because for 
vacuum oscillations y4£)_N oc RdP/dR and 6T oc i?c/P/(ii?, where the survival probability, 
P, depends upon the ratio of distance, R, to energy, E, i.e., P = P{R/E) [see discussion 
following Eq. (|A23|) and Eq. ( |A24| )]. It will be very difficult to test experimentally whether 



the relation given in Eq. (|32|) is present. 
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FIG. 9. The allowed regions for the [NC]/[CC] double ratio versus the day-night asymmetry, 
An-D) for a recoil electron energy threshold of 8 MeV. The meaning of the symbols is the same as 
in Fig. ^, except that the allowed regions now refer to [NC]/[CC] and tIn-d and a recoil energy 
threshold of 8 MeV. 

VI. [NC]/[CC] VERSUS DAY-NIGHT ASYMMETRY 

Figures |^ and ^ show maps of the currently- allowed regions in the Am^ — sin^ 26 plane 
onto the [NC]/[CC] — An-d plane for the electron energy thresholds of 5 MeV and 8 MeV, 
respectively. In Fig. ^ we show a simulated data point near the best-fit value for the LMA 
solution. The estimated la error bar for the [NC]/[CC] measurement is 3.6% after one year 
(Table II of Ref. 0) and is dominated by the statistical error in the determination of the 
neutral-current rate. For An-d, we assume for purposes of illustration a ±0.03 uncertainty 
in the absolute value as a la error, comparable to the accuracy that has been achieved after 
three years with the SuperKamiokande detector p^ . 
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The precision with which both the [NC]/[CC] ratio and the day-night asymmetry are 
measured will improve with time as more events are detected. 

Figure § and Fig. |^ are similar to Fig. ^ and Fig. because of the relation given in 
Eq. ( p2D between [NC]/[CC] and [CC]. The form of the allowed regions, the shape of the 
zones of avoidance, and the degree of overlap between different solutions are all quite similar 
in both sets of figures. But, because the double ratio [NC]/[CC] can potentially be measured 
with much better accuracy than [CC] alone, the correlations of [NC]/[CC] with An-d and 
other observables will have much stronger discriminatory power. 

A. Discriminating among solutions 

There are regions in the [NC]/[CC] — An-d plane in which only one oscillation solution 
is predicted to exist and, on the other hand, regions in which there is significant overlap 
between different oscillation solutions. There are also significant regions, easily visible to 
the eye in Fig. ^ and Fig. |^, in which no oscillation solutions are predicted to lie. 

We begin with a discussion of the regions where the identification of the oscillation 
solution may be unique and then discuss the ambiguous regions and the excluded regions. 

If An-d is observed to be greater than 0.2 and [NC]/[CC] is larger than 2.5, then the 
LMA solution will be uniquely singled out. The measurement of the first moment of the 
electron recoil energy distribution will provide a check on this inference since the LMA 
solution implies \6T\ < 0.01 (see Fig. ^ and Fig. |^), i.e., a very small distortion of the 
charged- current energy spectrum. Moreover, the reduced [CC] rate should be consistent 
with a value in the range 0.3 to 0.4 (see Fig. |] and Fig. ^). 

If [NC]/[CC] is measured to be larger than 4.5, then the only candidate solutions will 
be LMA and VACs. The two possibilities can be distinguished since (see Fig. |] and Fig. ^ 
LMA predicts a significant day-night asymmetry, An-d > 0.06, for large [NC]/[CC] and 
VACs predicts a very small asymmetry, |v4n-d| < 0.01. 

If, on the other hand, [NC]/[CC] is found to be smaller than 2.0, then the LMA and LOW 
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solutions will be eliminated. Values of [NC]/[CC] in the range 2.0 to 1.2 can be obtained 
with the SMA and vacuum solutions, but a value of [NC]/[CC] consistent with unity (and 
a small measurement error) would uniquely favor the MSW Sterile solution. In all cases of 
[NC]/[CC] less than 2.0, the predicted day-night asymmetry is very small, |An-d| < 0.02 
(see Fig. | and Fig. g). 

The most ambiguous region will be 2 < [NC]/[CC] < 4 and a small (< 0.02) day-night 
asymmetry. Figure |^ shows that multiple oscillation solutions can give rise to observables in 
this region. Moderate values of An-d (e.g., 0.02 to 0.12) and moderate values of [NC]/[CC] 
(e.g., 2.5 to 4.0) will also be ambiguous since all three of the MSW active solutions, LMA, 
SMA, and LOW can populate this region in the [NC]/[CC] — An_d plane. As we have 



discussed in Sec. |ITT|, a detailed study of the zenith angle distribution of the charged current 
events during the night may discriminate among these solutions. 

The zones of avoidance in the [NC]/[CC] — An-d plane are: all values of [NC]/[CC] larger 
than 5.2 (for any values of v4n_d, [CC], and 6T): An_d less than —0.02 (for any values of 
[NC]/[CC], [CC], and (5T); and [NC]/[CC] less than 2.5 together with An-d larger than 0.02. 

B. Correlation phenomenology 

The correlations between [NC]/[CC] and A^-d are similar to the correlations between 



[CC] and An_d that were discussed in Sec. |T|. The pointwise relation between [NC]/[CC] 



and ^N-D the LMA and LOW active solutions is washed out in Fig. || and Fig. ^ by the 
fact that v4^^j-) depends primarily on Am^ and [NC]/[CC] primarily depends upon sin^ 26*. 

In all cases, the day-night effect is small for vacuum oscillations. However, one can 
understand the general trend in Fig. and Fig. |^, according to which the solutions for 
VACl (larger Am?) correspond to smaller values of [NC]/[CC] than the solutions for VACs. 

It is easy to show from Eq. (g^), Eq- (^2^), and Eq. ( |A2^ ) that for the VACs solutions 



S = l + fc-AN-D, (33) 
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where k = /c(Am^) is a function of Am^ only. For the VACs scenario, there are four islands 
of solutions along which Am^ ^ constant but sin^ 29 changes. So, for a given VACs island 
k ^ constant, where k takes on a different value for each island of solutions. From Eq. (p3[), 
we see that [NC]/[CC] increases linearly with An-d- For An-d = 0, we obtain [NC]/[CC] 
= 1. All of these features are apparent in Fig. |^ and Fig. ^. 

The discussion following Eq. (|^) explains why for sterile neutrinos [NC]/[CC] is predicted 
to be close to, but not identical to, one. 



VII. [NC]/[CC] VERSUS SHIFT OF FIRST MOMENT 

Figures |10] and |ll| show, for electron energy thresholds of 5 MeV and 8 MeV, maps onto 
the [NC]/[CC] — 6T plane of the currently-allowed regions in the Am^ — sm^ 29 plane. For 



illustrative purposes. Fig. [1^ shows a simulated experimental point near the current best-fit 
predicted point for the LMA solution. The error bars are estimated la uncertainties from 
Table II of Ref. 0], with a 1.3% fractional uncertainty in the first moment and a 3.6% for 
the neutral to charged-current double ratio. 

Figure |10| and Fig. |ll| are similar to Fig. ^ and Fig. the latter pair refers to the 



correlation between the [CC] and the 6T variables. Because of the relation Eq. (|2^) between 
[NC]/[CC] and [CC], the acceptable solution space, as well as the zones of avoidance and 
the degree of overlap of the solution regions, are similar in the two sets of figures. 



A. Discriminating among solutions 

There are certain regions of the parameter space, [NC]/[CC] — 6T, that are relatively 
discriminatory. For example, for a 5 MeV threshold (see Fig. [1^), in the region in which 
6T > 0.03 and 1 < [NC]/[CC] < 5, only the VACs and the SMA solutions are represented. 
The most extreme values of the [NC]/[CC] parameter, e.g., [NC]/[CC] > 4.7 or < 1, would 
indicate, respectively, the LMA solution or the MSW Sterile solution. Only VACs solutions 
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FIG. 10. The allowed regions for the the double ratio, [NC]/[CC], versus the shift of the first 
moment, 6T, for a recoil electron energy threshold of 8 MeV. The meaning of the symbols is the 
same as in Fig. ^, except that the regions now refer to [NC]/[CC] and 6T. 

have a fractional shift 6T > 4.5% . Figure |T^ shows significant zones of avoidance, e.g., 
[NC]/[CC] > 2.2 and 0.01 < 6T < 0.01[NC]/[CC]. 

For an electron recoil energy threshold of 8 MeV (see Fig. ^), the degree of the overlap 
between the predictions of different scenarios is not larger than it was for a 5 MeV threshold, 
as is the case for [CC]—6T (cf. Fig. ^ and Fig. H). The allowed regions of VACs solution with 
largest Am^ do not overlap with LMA and LOW regions, as occur for a 5 MeV threshold. 
However, the largest Am^ VACs solutions do overlap with the allowed SMA region for an 8 
MeV threshold. 
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FIG. 11. The allowed regions for the the double ratio, [NC]/[CC], versus the shift of the first 
moment, 6T, for a recoil electron energy threshold of 8 MeV. The meaning of the symbols is the 
same as in Fig. ^, except that the allowed regions refer to [NC]/[CC] and 6T and the threshold for 
the recoil electron energy is 8 MeV. 

B. Correlation phenomenology 

For the SMA solution, both [NC]/[CC] and 6T are determined by a unique combination 
of neutrino variables, C,, defined by Eq. (PBD, so that the two measurables are strongly 
correlated: hst-nc/cc{^^-^) ^ 1- As the shift 6T increases, the double ratio [NC]/[CC] 
also increases (cf. Fig. |10| and Fig. pT| ). Using Eq. (^) and Eq. (|A12|) for the survival 



probability and Eq. ( [A16| ) for the shift of the first moment, we find for an electron threshold 
energy of 5 MeV 



[NC] 
[CC] 



^ e^^^ , (34) 
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where B = 29.3. 

For the VACs solution, a strong correlation exists. Using Eq. (^) and Eqs. ( [A23[ ) and 
(^), we find 

+ (35) 

where k' = k'{Am?) is the function of the Am^ only. For a given VACs island, k' can be 
considered as a constant. Thus the double ratio [NC]/[CC] is proportional to the shift 5T. 
The slope k' is different for different VACs islands. 

For the MSW Sterile neutrino solution, the rate [NC]/[CC] is strongly restricted by 
the measured value of [ESjgp^, whereas 5T varies over a significant range independent of 
[NC]/[CC]. 

The LMA solution does not predict a strong correlation, as can easily be seen from 
Eqs. (|A8|) and (^). The ratio [NC]/[CC] depends strongly on sin^ 26*, whereas 5T depends 
strongly on Aim?. The situation is similar for the LOW solution. 

VIII. CONCLUSIONS 

We discuss and comment on in this section the principal results from our analysis. We 
begin in Sec. |V111 A| with a restatement of the problem we address and then describe in 



Sec. [VIII B| our most important numerical results. In Sec. [VIII Q , we summarize how the 



predicted correlations and zones of avoidance between neutrino measurables can enhance 



the discriminatory power of solar neutrino experiments. Finally, we discuss in Sec. [VIII D 
some additional work that needs to be done in order to identify the correct set of mixing 
angles and mass squared differences of solar neutrinos. 

A. Correlated or not? 



For a given pair of neutrino parameters, Am^ and sin^ 26', the predictions for all the 
neutrino observables (e.g., day-night asymmetry or charged current rate) are completely 
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determined. Thus on a point-by-point basis the predictions for all the neutrino measurables 
are fully correlated. But, the currently allowed oscillation solutions constitute islands of 
finite size in the space of neutrino parameters. 

The practical question one wants to answer is: For a specified set of allowed solutions 
(e.g., LMA or VACs), how well correlated are the predictions for different neutrino mea- 
surables? In other words, if one considers for example the predictions that correspond to 
all the allowed values for Am^ and sin^ 20 currently included in the LMA solution, will the 
predicted values of observables like the day-night asymmetry and the charged current rate 
be strongly correlated? Or, will the range of Am^ and sin^ 26 within the allowed LMA 
domain obscure the point-by-point correlations? 



B. The answer 



Figures ^ to |TT] show the extent of the predicted correlations between different neutrino 
observables in the SNO experiment. These figures present our principal quantitative results. 
The specifics of these correlations depend upon the data set used (which will evolve with 
time) and the specified confidence level (99% in this paper). 

We have considered the following pairs of neutrino measurables: An_d versus [CC] (day- 
night asymmetry for the charged current, charged current rate), 5T versus [CC] (shift of 
the first moment of the charged current electron recoil energy spectrum, the charged cur- 
rent rate), An-d versus 5T (day-night asymmetry, shift of first moment), [NC]/[CC] versus 
An-d (double ratio of neutral current to charged current rate, day-night asymmetry), and 
[NC]/[CC] versus 5T (double ratio versus shift of first moment). For each pair of neutrino 
measurables, results are given for two different electron recoil energy thresholds, 5 MeV and 
8 MeV. The correlations are discussed in the text following the figures related to each pair 
of neutrino measurables. 

Some of the currently-favored neutrino oscillation solutions predict strong correlations 
among measurable quantities. For example, the allowed set of SMA solutions predicts strong 
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correlations between the values of An-d and [CC], as well as between ST and [CC] and 
between An_d and 5T . On the other hand, the LMA solutions predict correlations only 
between An-d and 5T and not between An-d and [CC] or between 5T and [CC]. 

The correlations, and the lack of correlations, can be understood from simple analytic 
arguments. We derive in Sec. [11 C| and in the Appendix approximate expressions giving 
the dependence of neutrino measurables upon /S.m? and sin^ 20. In subsections labeled 
'correlation phenomenology,' we describe the physical bases for the correlations and for the 
lack of correlations between different pairs of neutrino observables. 



C. Diagnostic power: correlations and zones of avoidance 

Does the simultaneous analysis of different observables enhance the diagnostic power of 
solar neutrino experiments? The answer is: "Yes, in some cases." In subsections labeled 
'Discriminating among solutions,' we emphasize for which cases the correlations among the 
predictions are strongest and how they can help in identifying the correct oscillation solution. 
We give examples in which multiple correlations can enhance the diagnostic power. For 
example, the values predicted by the LMA oscillation solution for the variables [NC]/[CC], 
An-D; and 5T are all correlated when An-d is large. We also show by examples that the 
dependence of the correlations on threshold provides additional constraints on the allowed 
solar neutrino solutions. 

The most powerful diagnostic pair that we have investigated may well be [NC]/[CC] and 
An-d- Figures |^ and ^ display the results for this case. This pair is particularly discrimina- 
tory because the systematic uncertainties in [NC]/[CC] and An-d can be reduced to values 
that are small compared to the ranges of the observables that are shown in Figs. ^ and |^. 
Moreover, correlations are predicted between the values of [NC]/[CC] and An-d for some 
favored oscillation solutions. By contrast, correlations involving the charged current rate, 
[CC], are severely compromised by the uncertainty in the value of the neutrino absorption 
cross section. Figures |TD| and |TT] also show significant correlations between [NC]/[CC] and 
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ST. 

Figures |^ to |Tl] show that there are zones of avoidance in the parameter space of neutrino 
measurables. None of the currently favored neutrino oscillation solutions predict values of 
the neutrino observables that lie within these unoccupied regions. We identify some of the 
more prominent zones of avoidance in the subsections 'Discriminating among solutions.' 

All of the currently favored oscillation solutions predict that the zones of avoidance will 
not be populated by values from experimental measurements. Thus an experimental test of 
whether or not the zones of avoidance are populated by future measurements is a general 
test of all of the presently allowed 2z/ oscillation solutions. 

D. Reducing the ambiguities 

Measurements with the SNO observatory will greatly reduce the allowed regions in neu- 
trino parameter space. Will a unique solution emerge from SNO measurements? Will we be 
able to identify the correct oscillation solution as one of the six currently-favored islands? 



There are many regions in Figs. to ITT] where multiple solutions (LMA, SMA, and 



LOW, e.g.) all overlap. In general, a unique identification will be possible only if one of the 
variables lies near an extreme value in one of the observables planes that we have considered 
in this paper. We will have to be somewhat lucky to be able to extract a unique solution 
from SNO measurements alone. 

In the future, we will study correlations between, on the one hand, measurable quantities 
in the SNO and SuperKamiokande experiments, and, on the other hand, quantities measured 



in low energy (less than 1 MeV) solar neutrino experiments (such as BOREXINO ||26|). We 
anticipate that unique inferences may be possible when low and high energy solar neutrino 
measurements are combined. 

Will the correlations and the zones of avoidance found in this paper also be valid for 
more complicated schemes of neutrino mixing three or even four neutrinos? Extensive and 
detailed computations are necessary in order to answer this question. 
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APPENDIX A: DEPENDENCE OF OBSERVABLES ON THE OSCILLATION 

PARAMETERS 

In what follows we present simple expressions for solar neutrino observables which are 
valid in narrow intervals near the best-fit points for different oscillation solutions. These 
expressions will be adequate for qualitative, and in many cases quantitative, understanding 
of the predicted correlations among the measurable quantities. Details of the approximations 



and more precise expressions for the observables are given elsewhere [E7 



We use results of numerical calculations of the iso-contours obtained in [15,16] in order 
to find parameters in the analytical expressions. The allowed regions in the Am^ — sin^ 29 
plane are taken from Ref. [jT^ (see Fig. |I|). 

Before proceeding to the approximate expressions valid for different oscillation solutions, 
we give the relevant definitions and equations that were used in deriving the results presented 
in the following subsections. 

For the MSW solution regions, the (daily) average survival probability is given by (see 



Eq. (35) in [|7 



P = ^{Pd + P^) = i [1 - cos2^5(l - 2P^)iUg - cos2^)] , (Al) 

where Pd and P/v are the averaged probabilities P(z/e — > i^e) during the day time and during 
the night time, respectively. 

The quantities that appear in Eq. ( |A1|) are defined as follows. 
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The probability Pi = P{ve ^i) is the probabihty that the solar neutrinos reach the 
surface of the Earth as the mass eigenstate Vi. 

The variable 6s is the matter mixing angle in the neutrino production region inside the 
Sun: 

(1 — 2ris cos 26* + ris) I 

where 

r/5 = ^m'/2EVs , (A3) 

and Vs is the matter potential in the center of the Sun; 
The regeneration factor, freg (see Eq. (30, 32) in 



describes the Earth matter effect. The quantity freg equals zero in absence of regeneration. 
In Eq. (P), 

riE = Am'^/2EVE, (A5) 

and Ve is the effective matter potential for the Earth. 

The Day-Night asymmetry is given by (see Eq. (37) in ||27|| ): 

_ Pn-Pd _ 2freg 

= - l/(l-2P0-cos2^ + /.e, ' ^^^^ 

where we have taken into account that in the Am^ region of significant Earth matter effect: 
?75 -C 1, and therefore cos205 ^ — 1. 



In calculating observables, the probability and asymmetry should be averaged over the 
neutrino energy. The effect of averaging can be represented by substituting for 2EV the 
effective parameters which are introduced in different equations below. The values of 
the effective m? should be determined from the results of exact numerical calculations. 
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1. LMA solution 



In the LMA solution region one has rjE ^ 1, so that according to ( [A4[ ) the regeneration 
factor can be approximated by 

sin^ 29 

Jreg ~ ^7" • (ATj 

Moreover, in this region Pi ^ cos^ 6 and rjs <^ 1. Then expanding cos2^5' given in Eq. 
( |A^ ) in powers of rjs and using the approximate expression Eq. ( |A^ ) for freg we obtain from 
Eq. (^H) the average survival probability 



P ^ sin^ e + - sin^ 26 
4 



A2 \ 2 2 
m \ m 



cos2e\ — ^ + 



E 



(A8) 



where m| and m'^ are fit parameters. The survival probability (and therefore the rate [CC] 
and the double ratio [NC]/[CC]) depend mainly on the mixing angle, sin^ 9; the dependence 
on Am^ is weak. The first term in the brackets of Eq. ( |A8[ ) is the correction due to effect 
of the adiabatic edge of the suppression pit, which is due to closeness of the resonance to 
the production point. This leads to deviation of a P from sin^ 9; m| is a parameter which 
corresponds to the product of the effective matter potential in the center of the sun, Vs, 
and the neutrino energy: m| ~ 2EVs. The second term in the brackets describes the earth 
regeneration effect, where m|; ~ 2EVe- 

From Eqs. (|l^) and ( |A8| ) , we find for the shift of the first moment 

sin^ 29 



5T 



2(1 - cos 2^) 



(A9) 



Here the negative term in the brackets is due to the adiabatic edge and the positive term 
describes the distortion due to the earth regeneration effect. The shift 6T is negative in the 
large- Am^ part of the LMA region and it becomes positive in the small Am^ part. For fixed 
Am^, the shift 6T increases with decreasing sin 29. For the reasons stated in the previous 
paragraph, 6T is very small for the LMA solution. 

For the day-night asymmetry, we find from Eqs. (^§) and ( [A7D the following analytical 
result: 
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A 



N-D 



Am? 



cos 29 m 
+ 



sin^ 29 



1 

2Am2 ^ 2 



(AlO) 



The asymmetry is, to a good approximation, inversely proportional to Am?. The first term 
in brackets leads to a decrease of the asymmetry when the mixing approaches maximal value. 
The second term is due to the regeneration effect in the earth. The last term describes the 
effect of the adiabatic edge which becomes important for Am^ ~ 10~^ eV^. For smaller 
Am^ the latter can be neglected and we obtain from Eq. ( AlO ) the equations for the iso- 
asymmetry lines: 



Am^ 



1 



A 



N-D 



sin^^ 
1 - cos 2^' 



(All) 



Comparing with results of numerical calculations, we find m| 6 x 10 ^ eV^ and 



ml ^ 3 X 10-6 eV^. 



2. SMA solution 



For not too small mixing angles, the survival probability can be well described by the 



Landau- Zenner formula [2S]: 



e ^0 . 



(A12) 



where ^ = Am^ ■ sin^ 29 and ^o(~ 2i?/ro) is a fit parameter, tq is the electron density scale 
height, 'n.e(r) oc exp(— r/ro). The effect of earth regeneration on the survival probability can 
be neglected here. From Eq. (|T3]) we find the reduced rate 

[ES]sK 



[CC] 



1 - r + re^/^o ' 



(A13) 



The day-night asymmetry can be parametrized in the following way: 

^1-2Plz(0' 



v4n-d = sin^ 29 ■ /(Am' 



az(0 



(AM) 



where 
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/(A™=)«.4-i^?;!4fi-, (A15) 

and A = 7.8 and ml = 3 x 10^^ eV^ are the fit parameters. Notice that ml corresponds to 
the Am^ with which neutrinos with an average detected energy resonate in matter of the 
earth. The probabihty -Plz(0 is given in Eq. (|A12|) . 



For the shift of the first moment, Eqs. (|T^) and ([A12|) yield 



6T = -Am^ ■ sin^ 29 = (A16) 
E E 



where is a fit parameter. 



3. LOW solution 

In the LOW region, the probability Pi equals the jump probability and can be approx- 
imated by the generalized Landau-Zenner probability. Pi ~ Plz- Also, in the LOW region 
t]e 1, and therefore the regeneration factor ( |A^) becomes 



/,,,^sin2 2^^ . (A17) 



Therefore the survival probability, Eq. (M), can be written using Eq. ( A17 ): 



P ^ sin^ 9 + - sin^ 29— —-r + P/ 7 cos 29 - ——-^ sin^ 29 . (A18) 



J ^2 



where {m'^) ~ 2EVe- Here the second term is the correction due to the earth regeneration 
effect which is important for the large Am2-part of the LOW region and the third term, 
which is proportional to the jump probabihty, Pl2 (defined below), is due to effect of the 
non-adiabatic edge. The generalized Landau-Zenner probability, P^g; which is valid for large 
vacuum mixing in the LOW region (see Ref. p9|), equals 



PLz ~ e"^"'"^ ^ , (A19) 

where 7 = 27TroAm'^ /2E and tq is the density scale height of the solar electron density 
distribution. 
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6T^ 2-— — -t + Plz^^ cos2^--— . (A20) 



From Eq. ([T7|) and Eq. ( [A18| ), we find for tlie first moment 

sin^2^ Am^ Am^ / ^„ Am^ 

4) "^L 

Here E/ttto) is the fit parameter. The first term in Eq. ( [A20| ) gives the effect of 

the non-adiabatic edge and the second term, which is proportional to Pl25 describes the 
distortion due to regeneration. For fixed sin^ 2^, the shift decreases with increasing Am^. 
In the small Am^ part of the allowed solution space, ST is positive since the spectrum is 
at the non-adiabatic edge of the suppression pit. The shift is zero at Am^ ^ 10~^ eV^ and 
then becomes negative due to the regeneration effect. The shift increases with decreasing 
sin^ 29. 

For the day-night asymmetry, we find 



Am 



1 - cos 29 Am2 2P/ 



+ 



LZ 



sm^29 2(m'^)2 (1 - 2Plz) sin^ 29 



(A21) 



The fit parameter {m'^Y = 2.5 x 10^^ eV^. In the region of the LOW solutions, the last 
term in brackets (containing P{^i) describes the effect of adiabaticity breaking and is positive, 
which suppresses the asymmetry. The asymmetry increases with Am^, in contrast with the 
behavior of the LMA solution. 



4. Vacuum oscillation solutions 

The standard expression for the vacuum oscillation probability leads to the following 
approximate relation: 

P = 1 - sin^ 29 ■ sin^ ^ . (A22) 

where my (~ AE/R) is a fit parameter. We find for the shift of the first moment: 

^T~pl^sin^2^.sinf^') . (A23) 
P m-Q \ my j 

The day-night asymmetry originates from the eccentricity of the earth's orbit and the 
existence of seasons [Q]. We find the residual asymmetry (after removal of the depen- 
dence of the total flux) which is related to the dependence of the oscillation probability on 
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distance from the sun. The stronger the dependence of P on the distance (oscillation phase) 
the larger the asymmetry. Clearly, the asymmetry is absent for P = constant. Therefore 



RdP 1 Am^ „ ^ Am^ 
— sin 26' -sin 2 — j- 



The expression in Eq. ( |A24| ) coincides with 6T, so An-d oc 6T. This result was obtained 



earlier (see the discussion following Eq. (^) using the fact that P = P{R/ E). 

For the VACl solution, there are strong averaging effects, and Eq. ( |A22D with a fixed 
characteristic energy does not reproduce accurately the functional dependence of the sur- 
vival probability on oscillation parameters. As a consequence, the relations Eq. ( A23 ) and 



Eq. ( |A24|) describe only very approximately the dependence upon neutrino parameters of 
the VACl solution. 

5. MSW Sterile solution 

The rate for the MSW Sterile solution is essentially fixed by the measured Su- 
perKamiokande rate, [ES]sk- The distortion of the electron recoil energy spectrum and 
the day-night asymmetry are similar to that for the SMA case, but the earth regeneration 
effect is much smaller. 
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